Stability properties of trapped Bose-Fermi gases mixture 
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The stability of Bose-Fermi gases trapped in an isotropic potentials at ultracold temperature is 
strongly influenced by the interaction between the fermions and the bosons. At zero temperature, 
the stability criterion is given in this paper using variation method, the results show that whether 
a fermion-boson mixture is stable depends mainly on the interaction between the fermions and the 
bosons. For finite temperature, however, the stability is not only related to the coupling constants, 
but also to the temperature. The stability conditions for finite temperature are also derived and 
discuss in details in this paper. 
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Since the realization of dilute alkali atomic vapor condensates(Bose-Einstein condensation or BEC) in 1995[1], large 
efforts have been make to study many-body effects and macroscopic properties of the gases, which may be more 
transparently demonstrated in BEC than in other many-body systems. For fermionic atomic vapor, however, it is 
difficult to achieve a degenerate gas. Since the evaporative cooling of a pure fermionic gas is ineffective at temperature 
sufficiently low due to the suppression of s— wave scattering between identical fermions. As theory and experiment 
advanced, a new rich phenomenology has appeared in which new conditions arise, which are not accessible in other 
BEC systems. One of the most stunning of these is the recent experimental demonstration of a condensate mixture 
composed of two spin states of a7 Rb[2]. The realization of two condensates mixture is related to the sympathetic 
cooling mechanism, i.e., the exchange of energy due to elastic collisions between atoms of cooled and thermal samples. 
Most recently, B.DeMarco and D.S.Jin [3] report their observation of degenerate Fermi gas using an evaporative 
cooling strategy. Although the strategy uses a two-component Fermi gas, the mixture of Bose and Fermi gas attracts 
a lot of attention from the viewpoint of both experiment and theoretical study. 

The mixed system of Bose and Fermi particles is itself an interesting subject for investigation. The hydrogen 
deuterium system has been studied at the early stage of these investigations [4], and there is now a lot of literature 
devoted to the properties of pure degenerate trapped atomic Fermi gases [5-9]. 

In a recent paper, Maimer has used a simple mean field models to study the spatial distribution of a Bose-Fermi gas 
mixture at T = OK within Thomas-Fermi approximation. The results show that the distributions depend strongly 
on the relative sign and magnitude of the boson-boson and boson-fermion scattering lengths. Here, we shall study 
the Bose-Fermi gas mixture using a variation method at zero temperature, this method was first introduced in[10] to 
study the BEC ground state in a harmonic trap of a Bose system, and later generalized by H.Shi and W.M.Zheng to 
study BEC with attractive interactions [11]. In addition, we study the stability of the Bose-Fermi gas mixture at finite 
temperature. The results show that there is a region of temperature in which the phase separation of the mixture 
happens. And the span of the region depends on the coupling constants. 

To begin, we consider a second-quantized grand canonical Hamiltonian of interacting Bose and Fermi gases 

H = H b + H f +V bf , 

Hb = J dr ^ + ( r )(^"^ + ^m b u; b r 2 )(t)(r) + ^f- J J drdr (f> + (r)(f) + (r )<f>(r )<j)(r), 
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H f = J drip + '(»")( _ M/ + 2 m / w / r2 )^( r )> 

Vbf = g b f J drdr <j) + (r)ip + (r )S(r - r )ip(r )4>(r), (1) 

where (j)(r) and ip(r) denote boson and fermion field operators with masses m& and m/, respectively. For weakly 
interacting dilute gases, the interactions between the bosonic atoms are modeled by S potentials and the interactions 
among the fermionic atoms are neglected, since the interactions between atoms at very low temperature is suppressed 
for polarized systems. g bb and g ba stand for boson-boson and boson-fcrmion coupling constant, respectively. 
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o-bb(o-bf) are s-wave scattering length between boson and boson (boson and fermion), and rribf is a reduced mass of 
the boson and the fermion. The chemical potentials Hb and [if are determined through the conditions 



N b = (J dr<j) + {r)(t>{r)),N } = { J dr^+ (r)V(r)). 



(2) 



At T = 0, self-consistent mean field theory, assuming that all N bosonic particles in a gas populated the same state 
denoted by single particle wave function $(r), lead to a nonlinear Schrodinger equation (or the Gross-Pitaevskii 
equation) for <£(r) = (<p(r)) 
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■ 9bbn b {r)]^{r) = E b ${r), 
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we here omit quantities gbfUf(r), which is smaller than gbbn-b(r) in the case of Nb >> Nf. In order to get a degenerate 
fermionic gas, the boson particles appear in the system only as a coolant, so the number of bosons is always much 
larger than the number of fcrmions. In the same approximation, the fermionic wave function is given by a Slater 
determinant 



*(ri,r 2 , ...,r Nf 



*i(ri) 
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where \Jij(r) is the single particle states determined by Hartree-Fock self-consistent equation 



The density of the fermions is given by 



fujjr 2 + g bf n b (r)]Vi(r) = E^i(r). 



n f (r) = |*(r)| 2 . 



(5) 
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In the semiclassical (Thomas-Fermi) approximation, the particle are assigned classical position and momenta, but the 
effects of quantum statistics are taken into account. Under this approximation, the Eqs.(3) and (5) for the boson and 
fermion wave function are equivalent to[13,14] 



-m b ^b r + 9bbn b (r) = fj, b , 



2m f 



[67r 2 n/(r)]3 + -mfujjr 2 + g bf n b (r) 



(7) 



The main conclusion of this equations is discussed in Ref.[13]. We obtain n b (r) = j^i^b — \'m b uj1r 2 ) from the first 
line of Eqs(7). Substituting n b {r) into the second line of Eqs(7), we yield 



^ r 9 / \ l — 99 tyb f / 1 9 9 \ 

[6tt n/(r)]3 + -m f uj f r + - — (p b - -m b uj b r ) = e F , 

2m f 2 g bb 2 



(8) 



this equation shows that the fermions experience a potential minimum in the center of the trap if g b f / g bb < 
mfLU 2 /m b uj 2 , in this case the entire distribution behaves like a fermionic core within the Bose condensate. The 
fermion density is a constant throughout the Bose condensate if g b f/g bb — mjUJ 2 /m b u 2 . Whereas the fermions are 
repelled from the center of the trap and localized near the edge: of the Bose condensate if g b f/g b b > m f^'f/ m b^ 2 - 
i.e. a phase separation occurs in this system. We would like to note that the distribution of BEC remains un- 
changed in the above discussions, since we assume N b » Nf. To drive Eq.(8), we assume that the Thomas- Fermi 
approximation(TFA) is valid. The coupling constant gbb and gbf may take any value as long as the TFA is available, 
and the phase separation depend mainly on ratio gbf/gbb- In what follows we discuss the separation of the bosonic 
and fermionic parts from the other aspect for zero temperature by using variation method, the results are indeed 
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different from those under TFA. We note the solution of Eq.(5) requires prior knowledge of the boson density profile 
n b = \$>(r)\ 2 . To obtain the density profile, we have to solve the Gross-Pitaevskii equation (3). There are a large 
number of literatures devoted to solve the Gross-Pitaevskii equation[12], we here use a variation method[ll] to solve 
the problem. For a isotropic trapping potential, we may assume the trial wave function for <!>(r) in Eq.(3) to be 

*(r) = ^N~ b Lui(^)h- m ^ r2/2h , (9) 
nn 

where u> is the effective frequency and is taken as a variational parameter. Substituting Eq.(9) into Eq.(3), we obtained 
the ground-state energy 

E b [$] = E b {u) = -N.hu + -N b h^- + g bb N 2 (^±)i. (10) 
4 4 uj Znn, 

If E b {(jj) is plotted as a function of uj , one sees that a stable local minimum exists only up to a certain maximum 
number of atoms for g bb < 0[11]. The critical point occurs where 

dE b (u). _ , d 2 E b {uj) 

— 0^ — \(u=u>c,N h =N bc ) = 0) and — — — \u=w c ,N h =N bc > 0. (11) 

Here, oj c stands for the variational parameter that minimizes the ground state energy. Using equation (10), for g bb < 
the critical number of bosons is given by 

7Vc <2 ^( Wc )-|J-(^)l. (12) 
1 9bb | m b 



where uj c satisfies 



hul - hJ b + 2g bb N b (^)iuj! = 0. (13) 



Parameter uj b ~ 166-ffz relevant to the experiment gives N b ~ 1400, which is in good agreement with the experi- 
ment^, 11]. The solution u> c of Eq.(13) against g bb is plotted in Fig.l, which shows that as \g bb \ increases, the variation 
parameter u c decreases, and it has a maximum equal to u b at g bb — 0. We will use this solution to study the stability 
of the mixture at zero temperature below. 

We may determine the ground-sate energy functional of the fermions provided n b (r) is known. In terms of the 
fermion distribution n/(r), the energy functional Ef of the fermions is given by [15] 

Ef = E f [n f (r)] = J ^ A-^n^r)] 5 / 3 + 1 J d 3 rm f u}r 2 n f (r) + J d 3 rg bf n f (r)n b (r), (14) 

since the interaction between the bosons and the fermions is rather week, we may consider a Gaussian function as a 
trial fermions' distribution 

n f (r) = ^^72(^)3/2^^^! _ (15) 
Here r/ and is treated as variation parameters. Substituting this wave functions into Eq.(14), one obtains 
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E f = E f (n,r f ,N f ,w c ) = P + ^nifUjjrjNf + hN b N f (irG)ie- Gr ~f (16) 
with 



P = P{Q,N f ) = Qi-(6TT 2 )zhnNf + 



4?7 



-Nf 



n{rrif\l + m b u c ) 

As known, a physical state corresponds to a stable or metastable point of the energy functional. If a separation of 
the fermion and boson component occurs, then rj c that minimizes the energy Ef takes a positive nonzero value. 
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In other words, there are no separations between the two components when energy Ef exhibits a minimal value at 
rf = 0. When distribution function is restricted to the form of the trial function (15) we may write the conditions of 
a minimal energy in terms of derivatives of the energy with respect to the adjustable variation parameters of the trial 
function. We show them as follows 

-n^ll -n 

an ln=n < ~ ' dr f Wf=rfc ~ ' 

d 2 E f d 2 E f . d 2 E f .„ „ 

( ana ) > °' ( 17 ) 



dfl 2 dr 2 K dfldr f 

the stationary conditions (there are no separation between the boson and fermion ) are 

Y = Y(N b ,u c ,g bf ) = lH-+ gbf N b m ^G(n c ,u; c ) — 2 + m ) ] 

x [m f uj 2 - 2g bf N b n 3 / 2 G 5 / 2 (n c ,iu c )} > 0, (18) 
where Sl c determined by c?i?/ / <9£l|r2 c = satisfies the following equation 

^(6^)1(3)1^1 _ l _ h q_ ? + gbfNh J ^n^ dG ^ <?£ = o, (i 9 ) 

the solution of Eq.(19) as a function of g b f is shown in Fig. 2, a magnification part of the curve near g b f = (but 
g b f > 0) is give in the inset. This curve indicates that the fermions prefer to occupy the trap centre for g b f < 0, and 
the larger the coupling constant 156/1(56/ < 0), the sharper the distribution of the fermions. However, as we show 
below, the fermions and the bosons can not always coexist even if g b j < 0. g^| r/=r/c = has two solutions, one 
solution is r/ c = and the another is 



r/ c = \J ^ln(2g b fN b iT2G2(Q c ,u) c )) - ln(m f ui 2 ). (20) 



For g b f < or 



mtui 2 f 

< 9bf < — 3^5 



2N b ^Gi{D. Cl uj c ) 

i.e. the interaction between fermion and the boson is attractive or weekly repulsive, the solution rf c = holds, which 
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indicate that there is not separation between the fermions and the bosons. For q b t > 3 m/ ^ f the fermions 

1 2N b n?G?(n c ,u> c ) 



experience a effective potential minimum at r/ c = y ^ln{2g b fN b 'KzGz (Cl c , u c )) — ln(mfU 2 ) > 0, the Bose condensate 

is surrounded by a shell of fermions in this case. Y(N b ,g b f,uj c ) as functions of the coupling constant g b f are shown 
in figure 3. We see that whether the boson- fermion mixture is stable depends not only on the coupling constant 
g b f and ^(through uj c ), but also on N b and ^(through O c ), i.e., the stability of the mixture system depends on 
the number of both boson and fermion system. For example, in Fig.3-a we show Y given by Eq.(8) as a function 
of the coupling constant g b f for fixed Nf — 100, N b = 10000, while Fig.3-b is for the same parameters as in Fig.3-a 
except for N b — 1000, it is obvious that the region of g b f in which the system has no phase separation has been 
broadened with N b decreases (for fixed Nf). The inset present the dependence of Y on g b f at a larger scale of g b f. It 
is interesting to compare the above mentioned results with those obtained by treating the fermions in the Thomas- 
Fermi approximation, this is done in Ref.[13,14], and we note that the semiclassical description gives a qualitatively 
correct description and it reliably predicts the phase separation. 

Now we tune our attention to discuss the above problem at finite temperature. First of all, we consider the 
homogenous case, for the boson and fermion system, thermodynamical properties are trivial if there are not interaction 
between them. But in this case the sympathetic cooling scheme does not take any effect and the degenerate fermions 
in a trapped potential have not been achieved. The thermodynamical properties may be changed when the interaction 
between the fermions and bosions is turn on, then a new phenomenon, the phase separation, may occur in this system, 
for a homogeneous fermion and boson mixture system, the Helmholts free energy can be written as[16] 

0F = + \9ffPfN f \ 2 f + ln(l - z b ) - ^ 9i {z b ) + 2g bbPb N b X 2 + g bf (X 2 + \})N f N b /V, (21) 
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where index / refers to the fermionic component, whereas index b stands for the bosonic one, iV, is the number of 
particles in component i, Xi denotes the thermal wave length of component i, f n (z) and g n (z) represent the Fermi 
and Bose integral, respectively. The equation (21) is based on the pseudopotential form of the atom-atom interaction, 
and may be assumed accurate when the system is dilute, i.e. pig\ << 1 and gu/Xi << 1, where pi is the density of 
the component i. This condition is well satisfies for the samples of alkali atoms in experiments to date[l, 12, 17,18]. 
From eq.(21) we obtain the chemical potential for each component straightforwardly, 



0p b = PfjP b + Ag bbPb X 2 b + g bf (X 2 b + X))N f /V, 
pp f = /Vf + 9ffPf>$ + 9bf{\l + \})N b /V, 



(22) 



where p® are the chemical potentials of ideal gas. There are three terms in each chemical potential, the second term 
comes from the interaction within the component and the third term is from the interaction between the fermion and 
boson component. As known, an homogenous binary mixture is stable only when the symmetric matrix p given by 



P : 



dpb d}J± 

dpb dp/ 

djjj_ djif_ 

_ dp b dp f 



(23) 



is non-negatively definite, in other words, all eigenvalues of matrix jj, given in Eq.(23) are non-negative. Mathemati- 
cally, for homogeneous fermion and boson mixture the stability conditions are 



^>0,^>0, 
dpb op f 



(24) 



and 



del 



dpb 9p b 

dp b dp/ 

dpj_ dpj_ 

dPb 9p f 



> 0. 



(25) 



For ideal gas, we have p b 



m{zb), Pf 



jrfz(zf), this leads to 



dp f fdz f y d Pb gdz b y 



(26) 



It follows from eqs (24) and (25) that 



Z(T,g bf ,g ff ,g bb ) = Z = (Ag bb Xl 



3//A/ 



9i{zb) 

h(zf) 



>o, 



and 



(27) 
(28) 

(29) 



It is well known that a homogeneous imperfect gas with attractive interaction is not stable. The fermions in this kind 
of gas could form BCS state, which consists two fermionic particles interacting with each other but not with the other 
fermions from the Fermi gas, whereas bosons with attractive interaction could collapse into liquid. Hence, we here 
discuss the system with repulsive interactions. It is obvious that the stability condition (27) and (28) hold always for 
g bb > 0, gff > 0. We would like to point out that the stability conditions (27-29) do not involve the densities of the 
both components. At first sight, this seems to be confusion, in fact, there is no contradiction. One can demonstrate 
that at low density the Helmholtz free energy of the bogoliubov gas reduce to a quadratic form in N b and Nf. To 
have a minimum, this form should be positive definite, i.e., det | d ^ b g Nf \ \ > 0. Therefore, the corresponding stability 
criterion involves only density-independent constants in the order of approximation used. This criterion is similar 
to the stability conditions for two-component Bose-Einstein condensate in a trapped untracold gas [19-25]. When 
T — > oo,Ai — > 0, hence Z ~ j^j- Thus at high temperature, the homogeneous binary gas mixture is always stable 
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and no phase separation occur. In the case considered here, Fermi temperature Tp — 2 mk B (f^y) ^ 1S mucn lower 
than BEC temperature T c = 2 J^ nKs ( 2 6iW ) ^ ' i- e -> as temperature decreases, it first passes the BEC transition point 
T c . When T T c , 51 (1) -» 00, so B 

A 3 . 

Z(T,g bf ,g bb ,g ff ) ~ 4 566 A^( 5// A^ + jT^y) - S6/(A? + A 2 ) 2 . 

In particular, when T « Tp, i.e., the temperature is much smaller than the Fermi temperature of the fermion system, 
the stability condition becomes (setting to/ = m b ) 

9bb9ff ~ ff6/ > 0, (30) 

which does not depend on temperature and coincides with the stability conditions of two-component BEC[19,25]. 
Although it is difficult to reach this region of very low temperature, yet it attracts much more attention. Because both 
superfluidity and shell effects are expected to occur at temperature much smaller than the Fermi temperature [6, 26]. Z 
given by Eq.(29) as a function of the temperature is shown in Fig. 4, we see that the system is always stable when T — > 
and T — > 00, and the system is unstable for T& < T < T c2 , where T c \ and T c2 are roots of Z(T,g bb ,g b f,gff) = 0. 
In particular, T c \ and T c2 depend on g b f g/f and g bb . As g b f decreases (for fixed g bb and 5//), T c \ tends to T c2 
(in Fig. 4 going from dotted line to solid line). The critical temperature T c \ and T c2 characterize the onset of the 
phase separation, which is quite different from the Bose-Einstein condensation and the degenerate fermions. The 
critical temperature of BEC and of the onset of degenerate fermionic gas depend mainly on the density of the system 
Ni/V(i = b, /). Especially, the BEC and the degenerate fermionic gas may happen even if g b f = 0. For the phase 
separation, however, nothing will happen if g b f = 0. For a fixed temperature and the coupling constant g b /,Z vs. g bb 
and g/f is shown in Fig. 5, which represents the dependence of the stability on the interaction strength inside each 
component. 

Until now, we considered only a homogeneous Fcrmi-Bose gas mixture at finite temperature. In reality, however, 
experiments with ultracold atoms are performed by trapping and cooling in an external potential that can be generally 
modeled by an isotropic harmonic oscillator V(r) = y^r 2 , where u>t is the trapping frequency. An exact criterion for 
the stability of an inhomogencous Bose-fermi mixture should involve calculating the Hclmholts free energy as a function 
at all eigenstates of the trapping potential. Fortunately, in the system considered here it is a good approximation to 
take use of the local-density approximation, which treats the system as being locally homogeneous. This requires that 
the level spacing hu>t of the trapping potential is much smaller than the Fermi energy. Of course, the local density 
approximation always breaks down at the edge of the gas cloud where the density vanishes and the effective Fermi 
energy becomes zero. In this approximation, the stability conditions can still be calculated by means of the equations 
derived above, with the understanding that now the effective chemical potentials are spatially dependent through 

M6 = M° - \m b ujlr 2 , n f = n° } - ^nifujjr 2 . 
Thus a local stability condition is the same as given in Eq.(29) but replacing Zi(i = 1,2) by 

2i = Zl e-^ mb ^ r \ and ~z 2 = z 2 e-* m ^>\ 

As shown in inset of Fig. 4, the region of temperature in which the system is unstable decrease for the case of r 7^ 0. 
As compared with the case without trapped potential, the total energy of the system increases for it in a trap. 
Alternatively, within the TFA, the chemical potentials decrease in this process. So this effect is equal to be that the 
particle number of the system has a loss. In this sense, the system is more stable than before. 

In summary, we considered a dilute Bose-Fermi gas mixture in an isotropic trap. The atom can interact via s-wave 
scattering except within the fermions. These interactions strongly affect the stability of the system at zero and finite 
temperature. In addition, the stability conditions depend on the ratio rate N b /Nf, the larger the ratio rate, the 
smaller the region of stability. For finite temperature, however, the stability conditions depends not only on the 
interactions, but also on the temperature. The region T c \ < T < T c2 in which the system is unstable depend on the 
strength of the interaction between and within the bosons and the fermions. For an anisotropic trap, the stability 
conditions remain unchanged, whereas somewhat would be changed for zero temperature compared with the case 
of isotropic trap. To study the effects, we should introduced the another variation parameter in Eqs (9) and (5) to 
characterize the BEC and the fermions in this trap. Consequently, the stability condition (18) for zero temperature 
changes and the phase separation could different for different orientation. These need further investigations. 
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Figure captions: 

Fig. l:The parameter ui c which minimizes the energy functional versus the coupling constant gbb- The trapped fre- 
quency ujb — 166-ffz and the number of the bosonic atom Nb = 1000. 

Fig. 2:The parameter f2 c which minimizes the energy functional as a function of the coupling constant gbf. The 
parameters chosen are gbb = 0.05 in units of huifa 3 (a = ^Jh/uibmb) and all the coupling constants are chosen in this 
units hence forth, Nf = 100, uif — 166-ffz. Scatter and solid line correspond to different number of bosonic atom, as 
specified in the figure. The inset presents the enlarged part of the curve near g^f = 0(> 0). 

Fig.3:Plot of Y given by Eq.(18) as a function of coupling constant gbf. The parameter chosen are &:Nf = 100, Nf, — 
1000. b:Nf — 100, Nb — 10000. The curve for a larger scale of gbf is presented as an inset in the figure. 

Fig.LPlot of Z given by Eq.(29) as a function of temperature T. The parameters chosen are Nf, = 1000, Nf = 
10000, 3 6f , = 0.05, 3// = 0.01. Dashed-dotted line:g 6/ = 0.3, dotted line g bf = 0.02, solid line g bf = 0.01. The dotted 
line in the inset is the same as the dotted line in the figure, while the solid line in the inset is for the gases in a trap 
with trapped frequency 166Hz. 
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Fig.5:Plot of Z as a function of and 56b. The parameters chosen are temperature T = O.lTp, gbf = 0.2. 
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